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U 

C^ Abstract 

'^N The response of bright and dark solitons to rapid variations in an ex- 

r-~^ pulsive longitudinal trap is investigated. We concentrate on the effect of 

^-H transient changes in the trap frequency in the form of temporal deha kicks 

and the hyperbohc cotangent functions. Exact expressions are obtained 

' , , ' for the soliton profiles. This is accomplished using the fact that a suitable 

^ linear Schrodinger stationary state solution in time can be effectively com- 

,.1^ bined with the solutions of non-linear Schrodinger equation, for obtaining 

f^ solutions of the Gross-Pitaevskii equation with time dependent scattering 

• length in a harmonic trap. Interestingly, there is rapid pulse amplification 

j^ in certain scenarios. 

S PACS nuniber(s) : 03.75-b, 05.45.Yv 

I 

r^ I Introduction 

o 

, ^, Bose-Einstein condensates (BECs) [1] - [3] can be confined to one and two 

dimensions in suitable traps [1] - [B], which has opened up the possibility for their 

^-H technological applications. BEC on a chip is a two dimensional configuration 

^ Hi whereas cigar shaped BECs refer to quasi one dimensional scenario [5] - 

iJ [10]. In the latter case, the transverse confinement is taken to be strong, with 

QQ the longitudinal trap frequency having comparatively much smaller value. In 

pg this condition, dark [TT] - [T3] and bright [13] - [TS] solitons have been observed 

■ ' in repulsive and attractive coupling regimes respectively. In the context of the 

(^ BEC, bright and dark solitons respectively correspond to density lumps and 

00 rarefactions compared to the background. The corresponding mean field Gross- 

^^ Pitaevskii (GP) equation being the familiar non-linear Schrodinger equation 

^ (NLSE) in quasi-one dimensions ^20^ , the observation of localized solitons has 

• '-j demonstrated the ubiquity of solitons in diverse branches of physics, which share 

rS common non- linear behavior. 

jrt Very interestingly, in the cigar shaped BEC, a number of parameters like the 

scattering length and trap frequencies can be changed in a controlled manner, 
leading to the possibility of coherent control of solitons. The former can have a 
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temporal variation due to Feshbach resonance, a scenario which has been exten- 
sively studied in the literature |HT] - [33] • In particular, the discovery of bright 
soliton and soliton trains has led to a systematic investigation of the impact of 
the change in the scattering length on the BEC profile. Recently, the temporal 
modulation of the non-linearity in one dimension arising due to the correspond- 
ing variations in the transverse trap frequency have been systematically studied 
in the context of generation of Faraday modes [21] , [2S] . 

In a recent paper [26], a method to obtain analytical solutions of the GP 
equation for a range of temporal variations of parameters like the scattering 
length, the longitudinal trap frequency etc. was demonstrated. It was shown 
that in the process of reducing the GP equation to a known non-linear equa- 
tion, having Jacobi elliptic functions as solutions, the consistency conditions 
obtained mapped onto the linear Schrodinger equation. This mapping enables 
us to obtain exact solutions of the GP equation for a wide range of temporal 
modulations of the control parameters. These solutions exhibit soliton trains, 
localized solitons (both bright and dark) for different regimes of the coupling 
parameter. Interestingly, exact solutions of the GP equation have been obtained 
for other kind of trap potentials, for example linear traps [2Z] using Darboux 
transformations. The effect of the variations in the time dependent coupling on 
the soliton profiles has also received a lot of attention [27] - [25] . 

The goal of the present paper is to investigate the soliton dynamics in a 
BEC in an expulsive harmonic trap, whose frequency is time dependent. The 
expulsive trap accelerates the soliton. Very interestingly it is found that for 
certain conditions it can be compressed, particularly under a sudden transient 
variation in the trap frequency. As we have mentioned earlier, a temporal mod- 
ulation in the transverse frequency leads to the generation of Faraday modes in 
the longitudinal direction. Hence, it is of interest to study the effect of transient 
variations in the longitudinal trap on the BEC profiles. We concentrate here 
on transient variations and for explicitness we use the delta functions and the 
hyperbolic cotangent functions to model rapid temporal variations in the trap 
frequency. We obtain exact solutions and study the spatio-temporal dynamics 
of the solitons and their response to the changes in the trap frequency. We point 
out that harmonic traps decorated with delta functions and a dimple potential 
have been recently studied numerically j30j and the properties like atomic den- 
sity, chemical potential, critical temperature etc. of the condensate have been 
obtained. 

In the following section, we obtain the exact solutions of the GP equation 
and illustrate how we can modulate the expulsive trap frequency using suitable 
quantum mechanical models [2S]. In section III, the trap frequency is tran- 
siently changed by hyperbolic cotangent functions. In section IV, we study the 
response of the solitons to sudden changes in the trap frequency in the form 
of delta kicks. For all the cases studied we obtain profiles of both dark and 
bright soliton trains and localized soliton solutions. The spatio-temporal dy- 
namics of the soliton solutions and the behavior of the non-linearity with time 
in the attractive and repulsive coupling regimes are plotted. It is found that 
for certain examples there is an amplification in the soliton amplitude resulting 



in its compression while in other examples the solitons spread with a decaying 
amplitude. 

II The Dynamics of the Bose-Einstein condensate 

For a BEC at zero temperature, confined in cylindrical harmonic trap Vo{x, y) 
m'w^(x^ + ?/^)/2 and time dependent harmonic confinement along the z direc- 
tion Vi(z,t) = m'uj},{t)z'^ /2, the dynamics is governed by the GP equation 

^^^-^ = {-^'^' + um^^t)\^+y)nv^t), (1) 

where V = Vb+^i, in' is the mass of the particle and U — 'i'K'h'^as{t)/m! gives the 
non- linear coupling with as{t) denoting the scattering length. The nonlinearity 
can be either attractive (a^ < 0) or repulsive (a^ > 0). In order to reduce the 
3-D GP equation to quasi-one dimensions, it is assumed that the interaction 
energy of atoms is much less than the kinetic energy in the transverse direction 
[S]. Substituting the following trial wave function in dimensionless units, 

1 , ( z \ ( . x^ + y'^ 



we obtain the quasi- ID NLSE 

idti^ = -\dzz4' + tWIV'I'V' + ^M(i)zV. (3) 

Here, 7(i) = 2as{t)/aB, M{t) — ajg(f)/w^, a_i = {h/m'uj±y^^ with as being 
the Bohr's radius. We have kept M{t), related to the trap frequency, time 
dependent. This enables us to modulate the trap frequency, allowing us to 
study the effect of different temporal profiles of the trapping potential on the 
solitons. The specific case of M{t) being a constant, describes the scenario 
where the oscillator potential is either confining (A/ > 0) or repulsive {M < 0). 
Analytic solutions to the above NLSE are obtained using the ansatz solution 

^iz,t) = ./A(r)F{A{t)[z-l{t)]}eMMz,t)h (4) 

where A{t) describes the amplitude and l{t) = L v{t')dt' determines the location 
of the center of mass. We assume the phase to be of the form 

$(z,i)=a(<)-ic(i)z2. (5) 

Substitution of V'(^,i) in the NLSE gives us the consistency conditions 

Ait) = ^0 exp ( I c{t')dt'\ , ^ + c{t)l{t) = (6) 

and 

7W=7o^, (7) 



where Ao,7o and Iq are constants. We also obtain a{t) = ag + ^5^ /^ A'^{t')dt' 
along with the Riccati equation for c(i), 

^-c\t)^M{t). (8) 

From the above equations, it can be seen that A{t),j{t) and l{t) are all non- 
trivially related to the phase component c{t). Added to this, the fact that the 
Riccati equation can be mapped onto the linear Schrodinger eigenvalue problem 
with potential V{t), 

- ^"{t) - M{t)(t){t) = (9) 

with M{t) — Mq — V{t), using the change of variable, 

c(i) = -|ln[0(t)] (10) 

gives us control over the dynamics of the BEC. Hence, knowing c(t) or the 
solution (f>{t) of ([9]), allows us to obtain analytical expressions of the control 
parameters through equations (10) and (|6|-(l7|. It is a well known fact that for 



many potential models, the Schrodinger equation can be solved exactly. In the 
next sections, we exploit this and the connection between c(t) and 4'{t) to study 
the dynamics of the BEC analytically, when the oscillator frequency undergoes 
a rapid change. 

The use of the above consistency conditions along with the ansatz solution 
in the NLSE gives us the following differential equation for F, in terms of the 
new variable T = A{t)[z - l{t)], 

F"{T) - \F{T) + 2kF^{T) = (11) 

with 1^ = —-^ and the differentiation is with respect to T. The general solution 
of this differential equation can be written in terms of the 12 Jacobi elliptic 
functions sn (T, m) , en (T, to) , dn (T, to) etc where to is the elliptic modulus 
[31, . The solitons will arise in the limit ?7i — > 1 , where dn (T, to) = en (T, to) -^ 
sech (T) and sn (T, to) -> tanh(T). 

Thus for 7o < 0, bright soliton trains of the form 



?A(2, t) = y/MT) cn (T/ro, to) exp[i$(z, t)] (12) 

exist, where Tq = —771^0/70 and A = (2m— l)/rQ (these coefficients are obtained 
using equations (|6]) -Q. In the limiting case to = 1, equation (12 1 corresponds 



to a bright soliton solution. It has been seen that for 70 > 0, there exist dark 
soliton trains of the form 



iP{z, t) = y/MT) sn{T/To, to) exp[i$(z, t)], (13) 

with Tg = itiAq/jq, X = — (TO,+ l)/rQ and in the limit m ^ 1, the above solution 
corresponds to a dark soliton. 



In the following sections, we obtain the profiles of bright and dark soliton 
trains and solitons in a BEC in an expulsive oscillator potential, with experi- 
mentally achievable temporal variations of the trap frequency. As seen in ([8|, 
the oscillator frequency M{t) depends on c{t). We first concentrate on the re- 
sponse of the soliton profiles to the variations in the trap frequency when c{t) 
is equal to the hyperbolic cotangent functions. 

Ill Dynamics of the BEC with c(t) = ± Bcoth (t) 

In this section, we look at the examples where the trap frequency of the 
confining trap is modulated using c{t) = ± -Bcoth (t) , B > 0. Using the consis- 
tency conditions in (|6| and ([t]) we obtain the analytical solutions of the NLSE 
as shown below. 

case (a) 

We take c{t) = Bcoth (i), B > in ^ and obtain 

M{t) = -B^ - B{B + l)cosech^(t). (14) 

Since Mq = —B^, we are looking at an expulsive oscillator scenario. We point 
out that V{x) = i?(i3 + l)cosech'^(j:) corresponds to the supersymmetric Rosen- 
Morse potential with superpotential W = Bcoth (a;) [35]. The function V{t) 
has a singularity at t = which causes a sudden change in the trap frequency. 
Substituting for c{t) in equations (l6l and (JTl one obtains 



A{t) = Ao 



o(t) = ao ^^ "^^^° (sinh(i) 4i) 


(15) 


sinh(i)^^ ,,^ ^sinh(f)^^ ,,^, 
sinh(to)j '^(*^ = ^°Uinh(to)j ''(^) = 


^ /sinh(fo)\'' 
\ sinh (t) / 

(16) 



which when substituted in ( 12 1 and ( 13 ), gives the following solutions describing 
the bright and dark soliton trains 

■4,{z,t) = V^osinh (t) en ( Aqsinh (t) ^^ ~ ocoscc ^ ^^ ^rn j exp(z^(z, t)), (17) 

iP{z,t) = v^Aosinh (i) sn ( Aosinh (t) ^^ ~ oco&cc ( >\ ^^\ exp(z^(z, £)) (18) 

respectively. The dynamics of these solutions are plotted in fig 1(a) and 1(b), 
where the matter wave density |'i/;(2,i)P is plotted for varying t and z. We can 
see that the soliton trains get amplified with time as we move away from the ori- 



gin. In the limit to -^ 1, from ( 17 1 and ( 18 1 we get the bright and dark solitons 



plotted in fig 1(c) and 1(d) respectively. These solitons are highly localized and 
are seen to be diverging from the t — Q line in the t — z plane. By changing 
the initial time to some t^, instead of t = in the consistency conditions and 
varying ^o, the location of the center of mass can be changed. The non-linearity 



7(i) is plotted against varying t, for the attractive (70 < 0) and the repulsive 
(70 > 0) coupling regimes in fig 1(e) and 1(f) respectively and we see that for 
a given 70, the non-linearity does not change sign implying we are either in an 
attractive or repulsive coupling regimes. By changing the constants B, Aq, 70, 
Iq and io the amplification, location of the solitons can be changed. 

Case (b) 

For c(i) = -Scoth (i), B > 0, M{t) = -B^ - B{B - l)cosech^(i). Here again 
we are looking at the expulsive oscillator scenario and the potential B{B — 
l)cosech^(a:) corresponds to the supersymmetric Poschl- Teller potential [52] , 
Moreover, this case is interesting because for B — 1, this potential corresponds 
to the free particle problem which has been studied in [2S] . Substituting c(i) in 
equations ^ and (|7|, we obtain 

Ait) = Ao ('^^Y , 7(t) = 70 (''^Y , m = ^0 ^ ^^"'^'^ ^ " 



\ sinh (<) / ' \ sinh (i) / ' \sinh(to) 

' (19) 

As in the above case, we obtain the profiles of both bright and dark soliton 
trains and localized solitons which are plotted in figure 2(a) - 2(f). Although 
the function V{t) has a singularity at i = 0, like in case (a), the response of the 
soliton solutions is entirely different. We see that the soliton trains have a very 
high amplitude around i = and they decay as we move away from the origin 
i.e. the soliton profile spreads with time. The amplitude and the decay time 
can be controlled by changing Aq. The solitons plotted in figures 2(c) and 2(d) 
are again highly localized, but in this case converge towards the i = line in 
the t— z plane. From figures 2(e) and 2(f), we see that though the non-linearity 
does not change sign for a given 70 , it does increase rapidly in magnitude around 
t — Q, unlike the previous case were 7(i) varied slowly with time. 

IV Response of BEC to delta kicks in the trap frequency 

In this example, we vary the trap frequency with an attractive delta potential 
located at the origin and see how the solitons respond to a sudden kick. Here, 

V{t) = -VQ5{t), yo>0. (20) 

and the solutions of ^ are 

Ut) = y|iexp(-^t) (t>0) 

^cxp(|^i) (t<0) (21) 

Ut) = J^eM^t) it>0) 



^jcxp{-^t) (t<0) (22) 



with Mq — — Y- Ii^ order to obtain (p2{t), we have made use of the fact that 
the second solution of a second order differential equation can be obtained from 
the first solution using the relation <j>2{t) = 4>i{t) J [dt'{(f>i{t'))~'^]. We can write 



c(t) in terms of (j)(t) using (10) and obtain the consistency conditions in terms 
of (j){t) as 



A{t) = Ao 



^(0) 



7(^) = 70 



0(0) 



l{t) = lo 



Ht) 



(23) 



>(i) ' ''' ">(t)' '' V(o)- 
Now we proceed to study the soliton profiles using the two solutions given in 



(21) and (22 1 separately. 



Case (a) 

First we consider the solutions in (21 1 and from equation (10) we obtain c{t) 



Vb/2 for t > and c{t) = -VqI2 for i < 0. From (23), we obtain 



y^ 



K), 



^0, 



A{t) = Aoexp(^ltl) , 7(t) = 7oexp(^|t|) , l{t) = ?oexp(-^|t|) (24) 



and also find that 



a(i) = flo + 



(A- 



2Vr 



^^^" [exp(Fo|i|) 



(25) 



Substituting A(t\^(€) and l(t) in (12), we obtain the solutions describing 



bright and dark soliton profiles, which are plotted in figure 3 along with the 
plots showing the behavior of the nonlinearity. From the figures we can see that 
the sudden temporal kick to the soliton profile results in rapid amplification of 
the soliton trains over time. Similar to the case (a) of the previous section, the 
solitons are highly localized and diverge from the origin. By varying Vo, we can 
create a soliton with the required amplitude. The non-linearity 7(i) for this 
case turns out to be an exponential function. Thus the scattering length a^ii) 
turns out to be an exponential in time |18j . It is known that the dynamics of 
a BEC in an expulsive parabolic potential are governed by an NLSE with an 
exponentially varying scattering length. The conditions in which these systems 
are one dimensional and the range of the parameter values for which the con- 
densate is stable are discussed in [16,. 

Case (b) 



Here, we repeat the above analysis with the solution given in (22). For this 



solution, the nonlinearity 7(i) changes sign when we go from the negative time 
domain to the positive time domain. Here, we study the soliton dynamics for 
i > and the corresponding expressions are 



c(t) 



2 



a(t) = ao 






[exp(-yoi) - 1] 



Vb Vo Vo 

A(t) = AoCxp( — -t) , 7(i) = 7oexp(-yi) , l{t) = /oexp(-yi). 



(26) 



(27) 



The soliton profiles are plotted in figure 4. We observe that the soliton train 
amplitude decays rapidly as t increases unlike the previous case. For smaller 
values of Vq, the trains exist for longer time. Thus the decay of the solutions 
can be controlled by varying Vq. 

In conclusion, we have analytically treated the effect of variations of the 
longitudinal trap frequency on the BEC profile. We have looked at the cases 
where the trap frequency is modulated using hyperbolic cotangent functions and 
delta kicks. It is observed that in the gain/loss less scenario, these variations in 
the trap frequency resulted in the soliton profiles either amplifying or spreading. 
We have also shown that we can manipulate the amplification and the location 
of the solitons by varying the control parameters. 

The fact that a stationary Schrodinger eigenvalue problem can be combined 
with solutions of non-linear Schrodinger equation to obtain the solutions of 
Gross-Pitaevskii equation with time dependent parameters, enabled us to ob- 
tain the exact dynamical behavior. As mentioned in the introduction, there are 
other methods of obtaining exact solutions of the GP equation. The method 
used in this paper is simpler and the fact that the Schrodinger equation can be 
solved exactly for many potential models, offers us a huge choice of temporal 
variations of the trap frequency from the class of exactly solvable models. We 
can choose the models which can be realized experimentally. Along with this 
wide choice, the other advantage is that we have total control over the soliton 
dynamics. It would be interesting to see how solitons evolve in space and time 
for other choices of variations in the trap frequency. Another interesting prob- 
lem will be the study of two soliton dynamics under similar variations in the 
trap frequency [53] ■ 
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(a) Bright soliton train, m = 0.5 
TO = 0.375 




(b) Dark soliton train, m = 0.5, tq 
0.375 




(c) Bright soliton, m = 1, to = 1 
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(e) ■y{t) versus t, 70 = —0.5 



(f) ■y{t) versus t, 70 = 0.5 



Figure 1: c(i) = coth{t) in the expulsive oscillator regime. Plots of the 
bright [(a),(c)] and dark [(b), (d)] soliton solutions, ii{z,t) = yjA{t)F{{A{t)[z~ 
?(i)])/To,m}exp[za(t) — ^c{t)z'^] with F = cn{T/To,m) and F = sn{T/To,m) 
respectively, are shown. The plots (a)-(d) show \ip{z, i)p versus t and z and the 
plots (e) and (f) show the variation of the nonlinearity "fit) with t. The values 
of the constant parameters are Aq = 0.5, Iq = 5, to = 0.9 and B = 2. 
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(a) Bright soliton train, m = 0.5, tq 
0.375 




(b) Dark soliton train, m = 0.5, tq = 0.375 





(c) Bright soliton, m = 1, tq = 1 



(d) Dark soliton, m = 1, tq = 1 
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(e) ■y(t) versus i, 70 = —0.5 
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(f) 7(t) versus t, 70 = 0.5 



Figure 2: c(i) = —coth{t) in the expulsive oscillator regime. Plots of the 
bright [(a),(c)] and dark [(b), (d)] soliton solutions, ■ip{z,t) = y/A{t)F{{A{t)[z- 
l{t)])/TQ,m}exp[ia{t) — ^c{t)z'^] with F ~ en {T / t^ , m) and F — sn{T/To,m) 
respectively, are shown. The plots (a)-(d) show |^(2;,i)p versus t and z and the 
plots (e) and (f ) show the variation of the nonlinearity ^{t) with t. The values 
of the constant parameters are Aq — 0.5, Iq = 5, to = 0.8812 and B ^ 2. 
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(a) Bright soliton train, m = 0.5, 
TO = 0.375 




(b) Dark soliton train, m = 0.5, tq 
0.375 




(c) Bright soliton, m = 1, ro = 1 
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(d) Dark soliton, m = 1, tq = 1 









-OS 


^^^^ — 


^^--^^ OS 




-5 


^N 




M')-'" 


\ 




-IS 


\ 




-» 


\ 




^35 


\ 




(e) 7(t) versus t, 70 = —0.5 



(f) 7(4) versus t, 70 = 0.5 



Figure 3: Attractive 5 potential in the expulsive oscillator regime, (j){t) = 
VW2cxp(-V"o|i|/2). Plots of the bright [(a),(c)] and dark [(b), (d)] soliton 
solutions, i/'(;j,i) = yj A{t)F{{A{t)[z - ?(t)])/To,m}exp[ia(i) - \c{t)z'^] with 
F = cn(T/To,m) and F = sn{T/To,m) respectively, are shown. The plots (a)- 
(d) show |-(/'(z,t)p versus t and z and the plots (e) and (f) show the variation of 
the nonlinearity 7(t) with t. The values of the constant parameters are Vq = 10, 
Ao = 0.5, If) == 5. 
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(a) Bright soliton train, m = 0.5, 
TO = 0.375 




(b) Dark soliton train, m = 0.5, tq 
0.375 





(c) Bright soliton, m = 1, rg = 1 



(d) Dark soliton, m = 1, tq = 1 
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(c) 7(t) versus t, 70 = —0.5 



(f) 7(4) versus t, 70 = 0.5 



Figure 4: Attractive (5 potential in the expulsive oscillator regime, (j>{t) = 
y/2/Voexp{Vot/2). Plots of the bright [(a),(c)] and dark [(b), (d)] soliton 
solutions, ■ip{z,t) = y/A{t)F{{A{t)[z - l{t)])/To,m}exp[ia{t) - \c{t)z'^] with 
F = cn(T/To,m) and F = sn{T/To,m) respectively, are shown. The plots (a)- 
(d) show |-0(z,t)p versus t and z and the plots (e) and (f) show the variation of 
the nonlinearity 7(t) with t. The values of the constant parameters are Vq = 10, 
Ao = 0.5, lo = 5. 
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